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We theoretically examine entanglement in fractional quantum hall states, explicitly taking into 
account and emphasizing the quasi-two-dimensional nature of experimental quantum Hall systems. 
In particular, we study the entanglement entropy and the entanglement spectrum as a function of 
the finite layer thickness d of the quasi-two-dimensional system for a number of filling fractions v 
in the lowest and the second Landau levels: u = 1/3, 7/3, 1/2, and 5/2. We observe that the 
entanglement measures are dependent on which Landau level the electrons fractionally occupy, and 
find that filling fractions 1/3 and 7/3, which are considered to be Laughlin states, weaken with d in 
the lowest Landau level (j^=1/3) and strengthen with d in the second Landau level {i'=7/3). For 
the enigmatic even-denominator i/ — 5/2 state, we find that entanglement in the ground state is 
consistent with that of the non-Abelian Moore- Read PfafEan state at an optimal thickness d. We also 
find that the single-layer v = 1/2 system is not a fractional quantum Hall state consistent with the 
experimental observation. In general, our theoretical findings based on entanglement considerations 
are completely consistent with the results based on wavefunction overlap calculations. 

PACS numbers: 73.43.-f, 71.10.Pm 



I. INTRODUCTION 



The discovery of the fractional quantum hall effect 
(FQHE) in 198^ has proven to be one of the most 
significant experimental findings in all of physics. The 
incompressible quantum fluid that manifests in a two- 
dimensional (2D) electron system in the presence of 
a strong perpendicular magnetic field at low tempera- 
tures can not be explained by the "conventional" Lan- 
dau theory of phase transitions^. Instead this unique 
phase is an example of an emergent topological state 
of matter, where the charged quasi-particle excitations 
are governed by anyonic rather than fermionic or bosonic 
statisticd^''^. The unique nature of FQHE becomes ob- 
vious with the realization that the phenomenon occurs 
only in the truly strongly interacting limit of vanishing 
(or extremely small) kinetic energy with the noninter- 
acting ground state having a macroscopically large de- 
generacy. Although the FQHE was discovered decades 
agcPm^ many questions still linger to this day. The most 
notable of these is the nature of the FQH state observed 
at electron filling fraction v = 5/2. First discovered in 
198/221, the ^ 5/2 state is (so far) the only excep- 
tion to the famous "odd denominator" rule for mono- 
layer FQH systems given by the Laughlin a nsat and 
the more general "composite fermion" theorjl^HU. (The 
odd denominator rule assumes FQHE in a monolayer sys- 
tem. Bilayer FQH systems, however, have several even 
denominator states that are well-understood. See, for 
example Refs. fT2Hl4|) . Currently, the leading theoretical 
candidate for the v = 5/2 state is the Moore- Read (MR) 
Pfaffian model state introduced by Moore and Read in 
1990^^1 Recently it has been noted that the particle-hole 
conju gate to the MR Pfaffian, the so-called anti-Pfaffian 
liSlill, is also a viable candidate for the FQHE at = 5/2. 



An exceptional feature of these ansatz is the prediction of 
an yonic qu asi-particles with non-Abelian braiding statis- 
ticslSIiSIiSI This prediction is attractive giv en th e current 
interest in topological quantum computin^^HUI]^ ^jjg 
true nature of the FQH v = 5/2 state is, by no means, 
a settled questiorP^HUl. Directly probing the topological 
nature of FQH states has proven to be both a theoretical 
and an experimental challenge, making definitive and un- 
controversial verification of the MR Pfaffian ansatz elu- 
sive. However, recent developments in the field of quan- 
tum information have shown that measures of entangle- 
ment are useful tools in examining the global features of 
many-body strongly correlated quantum state^^. 

The true nature of the v = 5/2 FQHE state is one of 
the most prominent open questions in condensed mat- 
ter physics, and is a primary motivation for the current 
study. Recent experimental studies have explored this 
mysterious state and give some weight to the Moore- 
Read theory that is believed to explain it. One of these 
studies is the recent experiment by Venkatachalam et. 

that measured the charge of localized excitations 
in the v = 5/2 state to be e/4 as predicted by the 
MR theorjEH. These results are consistent with previ- 
ous studies by Radu et. al. and Dolev et. al. that used 
shot noise to investigate the local charg^^ZHMl, Jn addi- 
tion to the experiments noted above, Willett et. alW^ 
have seen evidence of quasiparticle interference oscilla- 
tions that support the existence of charge e/4 excitations 
at v = 5/2. Another recent study is work performed by 
Bid et. alW^ that experimentally observed the theorized 
neutral m ode o f the v = 5/2 state consistent with the 
MR theorjl221MI Although these developments point to 
the MR theory as the likely candidate for the v = 5/2 
state, they are not sufficient to unambiguously establish 
the existence of non-Abelian anyons. Also, it should be 
noted that the MR theory predicts a spin polarized state 
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at V = 5/2, but recent experimental worlP^l^fil suggest 
that this state may be unpolarized in some cases. Al- 
though we do not consider spin in this work, we attempt 
to leverage recent developments in quantum information 
theory to observe topological features predicted by the 
Laughlin and MR theory in numerically obtained FQHE 
ground states and examine how these features change un- 
der a realistic change (in particular, the finite thickness 
effect with varying thickness) in the effective interaction. 

These developments in quantum information theory 
provide the context for our current work. In this arti- 
cle we provide a detailed numerical theoretical study of 
entanglement for FQHE states incorporating the depen- 
dence on the quasi-2D layer thickness of the transverse 
dimension (i.e the finite layer thickness effeclP^^^. We 
emphasize physics of the finite-thickness effect in this ar- 
ticle partly because it provides a qualitative understand- 
ing of the FQHE in higher Landau levels as shown in 
Ref. riband Indeed, the orbital Landau level (LL) de- 
pendence of the FQHE is not completely understood'^^. 
The theory behind the MR Pfaffian model state, for ex- 
ample, makes no distinction between the half-filled sec- 
ond orbital Landau level (SLL) (i.e., ly = 5/2) and the 
half- filled lowest orbital Landau level (LLL) {i' = 1/2), 
but no FQHE has so far been observed at = 1/2 in 
monolayer systems. Furthermore FQH states in the SLL 
are relatively rare and generally "weaker" (i.e., requir- 
ing lower temperatures and higher sample mobilities to 
experimentally observe due to their relatively small en- 
ergy gaps on the orde r of 0.5 K or less) compared to 
FQH states in the LLlPHHESl. Some light was shed on 
this phenomenon of fragility of the FQHE in the SLL in 
Ref . l42l and l43l where it is shown that the finite-thickness 
effect is qualitatively dependent on Landau level. In par- 
ticular, it is shown that a finite, non-zero, layer thickness, 
d, helps stabilize the FQHE in the SLL, whereas in the 
LLL, finite d tends to "weaken" the state. Also, for the 
ly = 5/2 state, there appears to be an optimal thickness 
where the MR Pfaffian is "strongest." Thus, there seems 
to be a close connection between finite thickness and the 
SLL FQHE, which we explore in this work by calculating 
the thickness dependent entanglement properties of the 
FQHE. The results given in Ref. HH and H3] are primar- 
ily based on calculations of the overlap between numeri- 
cally obtained ground states and model FQHE wavefunc- 
tions^^ — in particular, the Laughlin and MR Pfaffian 
model states. Although the overlap is a powerful theo- 
retical tool — in fact, the wide acceptance of the Laughlin 
wavefunction as the appropriate description of the exper- 
imentally observed odd-denominator FQHE is arguably 
based on overlap results — it is not always definitive and 
can be misleading in some cases. An example of this in- 
volves the ly — 2/5 state that has been shown to have 
a large overlap with both the Jain composite fermion 
wavefunction as well as the so-called "Gaffnian" wave- 
function, even though these two states have different un- 
derlying topological order belonging to different univer- 
sality classes^^ . Therefore, we seek alternative theo- 



retical tools to probe the Landau level dependence of the 
FQHE through the finite-thickness effect, which should 
inevitably give insight into the nature of the enigmatic 
V ~ 5/2 state. In this article, we examine bipartite en- 
tanglement as an alternative measure (in particular the 
entanglement entropy and the entanglement spectrum) 
and study its dependence on finite layer thickness, d, 
in quasi-2D FQH states. One reason for our study of 
the thickness-dependent FQHE entanglement is that the 
thickness parameter d enables a continuous tuning of the 
Hamiltonian, lending to a continuous variation in the en- 
tanglement allowing a comprehensive systematic study. 

Bipartite entanglement measures are tools designed to 
quantify the extent to which degrees of freedom are en- 
tangled in a bipartitioned systerrP^. The most straight 
forward of these is the entanglement entropy (EE), de- 
fined as the the von-Neuman entropy 

S'b = Tr[/9AlnpA] (1) 

of the reduced density matrix pAjB = Tr5/^[|^') ('I'|] for 
state I vt) & H = Ha in a Fock space % that has 

been partitioned into two parts. The EE has proven to be 
a very powerful tool in examining quantum correlations 
in interacting many-body systemJ^S!. In particular, the 
scaling of the entanglement entropy with system size has 
been shown to follow certain "area" laws that can identify 
quantum phase transitions in certain cases^^. Also, the 
EE can be used to extract the "topological entanglement 
entropy" which is an indicator of topological order in 
the systerrPSiSZ' (i.e. states in the same topological class 
will have the same topological entanglement entropy). 
Given these developments, the EE and the topological 
EE appear to be attractive tools to probe FQH states. 
However, as we discuss below in more detail, obtaining a 
precise quantitative estimate of the topological entropy 
in FQH states requires a technical procedure prone to 
introducing significant errorJ^^^S^. In our study, we do 
not attempt to make such an estimate. Instead we focus 
on FQHE states with the primary objective of observing 
the qualitative behavior of the EE as a function of finite 
layer thickness, d, in the LLL and SLL. Our goal is to 
obtain the qualitative dependence of the EE in FQHE as 
a function of layer thickness d. 

We also investigate the "entanglement spectrum" of 
quasi-2D FQH states in this study. Introduced by Li 
and Haldane^^, the entanglement spectrum (ES) provides 
more information than the entanglement entropy alone. 
In particular, Li and Haldane conjectured that there is 
a direct correspondence between the low-lying eigenval- 
ues of the operator h = — ln[pyi] and the edge modes 
of the system and, thus, ES can be used to determine 
characteristics of the underlying conformal field theory 
(CFT) of the corresponding FQH ground state. As long 
as these "CFT-like" states are well-separated from the 
"generic non-CFT-like" states in the entanglement spec- 
trum of the ground state in the thermodynamic limit, 
then it is conjectured that the identified CFT does, in- 
deed, describe the state. In other words, states described 
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the by the same CFT (i.e. in the same universaUty class) 
will have the same low-lying structure in their respective 
entanglement spectra. Numerical studie^^^HMI and some 
recent analytical resultJ^ support the Li and Haldane 
conjecture. In our study, we carefully examine the entan- 
glement spectrum of finite sized FQHE states by quan- 
tifying the separations between the suspected low-lying 
CFT and non-CFT-like states ("entanglement gaps") as 
a function of finite-layer thickness. The entanglement 
gaps serve as a semi-quantitative measure of how well 
the state in a finite system fits with the universality class 
described by the CFT. By doing this we discover that 
the entanglement gaps follow trends qualitatively similar 
to the EE, leading to similar conclusions obtained in the 
overlap calculations given in Ref. [321 and 33] (i.e., that 
the finite-thickness effect strengthens the FQH states in 
the SLL). We also examine the entanglement spectrum 
in the so-called "conformal limif'^S. The aim of the con- 
formal limit is to remove finite size artifacts from the 
geometry of the system (in our case, the sphere), allow- 
ing for a concrete definition of a full entanglement gap. 
We find that the conformal limit does result in a full 
entanglement gap in most cases studied, but not in all 
situations. The lack of a entanglement gap in these ex- 
ceptional cases, however, is likely due to our choice of pla- 
nar pseudo-potentials in our worl^SS* further work along 
this line would be necessary to fully understand these 
situations. 

The structure of this article is as follows: In section 
\n\ we describe our methods for numerically obtaining 
the exact ground state of the Coulombic Hamiltonian 
for FQH ground states at filling fractions v = 1/3, 7/3, 
1/2, and 5/2 and our chosen model of the finite thickness 
of the quasi-2D system as well as the model wavefunc- 
tions (Laughlin and MR Pfafhan) to which we compare 
the numerically obtained ground state. Further, we de- 
scribe how our entanglement measures are defined and 
calculated. In section \IU\ we provide our results for the 
entanglement entropy (III A I, entanglement gaps in the 
entanglement spectrum ( III B ) , entanglement spectra in 
the conformal limit ( III C I and density of states calcu- 



lations of entanglement spectra ( III D I for FQH ground 
states at the filling fractions listed above. We provide 



our conclusions in section |IV| Finally in the appendix, 
we discuss our choice of planar Haldane pseudopotentials 
over spherical pseudopotentials and examine the implica- 
tions of this choice by comparing entanglement spectra 
in FQHE ground states at zero thickness obtained with 
either choice of pseudopotentials. 



II. METHOD 

We begin by considering a quasi-2D geometry where 
spinless electrons are confined in the x-y plane with layer 
thickness d along the z-axis and an external magnetic 
field B also along the z-axis. In the non-interacting case, 
the presence of the magnetic field quantizes the electron 



energy levels into highly degenerate Landau levels (LL) 
with energies En — {n + 1/2) TiWc where lUc = eB /mc and 
n = 0,1,2,... is a non- negative integer that defines the 
LL index. The degeneracy of each LL per unit area is 
given by (27tP)~^, where I = \fficfeB is the magnetic 
length and defines a length scale for the problem. The 
filling factor is given by = Ns/{'2,ttP) where is the 
particle density per unit area. When we include spin, 
the degeneracy of each Landau level is doubled. If we 
assume the electrons are polarized by the magnetic field 
and the LL's are filled sequentially by spin up and spin 
down electrons, then Q < v < 2 corresponds to states 
in the lowest Landau level (LLL) with LL index n — Q 
and 2 < J/ < 4 corresponds to the second Landau level 
(SLL) with LL index n = 1. Thus, the 1/2-filled LLL 
corresponds to v = 1/2 or 3/2 and the 1/2-filled SLL 
corresponds to v = 5/2 or 7/2. The identifications of the 
1 /3- filled filling factors is obtained similarly. 

The presence of electron-electron interactions clearly 
complicates this picture, but if we assume the electrons 
are confined to a single LL (i.e., there is no LL mixing), 
then the kinetic energy is a constant that can be ignored, 
giving the effective Hamiltonian: 



N 



H 



(2) 



where N is the total number of particles and Tjj — \ri^rj\ 
is the distance between particles measured in units of the 
magnetic length Note that our assumption of no LL 
mixing might not be a very good approximatiorP^HHI \y\ 
all cases especially when considering the FQH state at 
V — b/2 where it has been observed at "low" magnetic 
fields where the LL mixing parameter [{e^ / {el))]/fujjc ^ 1 
{e^ jd characterizes the strength of the Coulomb interac- 
tion where e is the electron charge and e is the dielectric 
constant of the host semiconductor). Furthermore, the 
Hamiltonian we are using is completely particle-hole sym- 
metric and apparently the ground state does not sponta- 
neously break particle-hole symmetry^^. This means that 
all of our conclusions about the MR Pfaffian state apply 
equally to the anti-Pfaffian state since the two ansatz 
are particle-hole conjugates. LL mixing, as shown re- 
centljCsHZH, 

can induce three-body terms which explicitly 
break particle-hole symmetry leading to a possible pref- 
erence toward the MR Pfaffian or anti-Pfaffian, however, 
in this work we do not consider such terms and can make 
no distinction. 

This Hamiltonian (Eq. ([2|) can be parametrized 
in terms of the relative angular momentum between 
two particles jxiij by the Haldane "pseudopotential" 
expansiorP^ 



E 

m{odd) 



N 



(m,j) (3) 



where Pm{fnij) projects onto states with relative angu- 

(n) 

lar momentum rriij = m and Vm are the Haldane pseu- 
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dopotentials for a given relative angular momentum m 

(n) 

and Landau level index n. These pseudopotentials Vm 
are the energies of a pair of particles with relative angu- 
lar momentum m confined to the n-th LL. For spinless 

(n) 

fermions, only Vm for odd m enter the Hamiltonian due 
to the Pauli exclusion principle. In a planar geometry, 
the Haldane pseudopotentials are given in terms of the 
Fourier transform of the interaction potential V{k) by 

Vi^^ = / dfcfc[L„(fcV2)]2L,„(fc2)e-feV(fc) (4) 
Jo 

where Ln{x) are Laguerre polynomials. This expansion 
allows us to work strictly in the Hilbert space of the low- 
est Landau level since all necessary information about 
electrons confined in higher Landau levels are contained 
in the pseudopotentials. 

The model wavefunctions (Laughlin and MR Pfaffian 
states) can be obtained by diagonalizing certain "hard- 
core" Hamiltonians. The Laughlin wavefunction 



N 



Laughlin 



i<j 



(5) 



where z = x—iy is the electron coordinate in the complex 
plane, and is the zero-energy ground state of a special 
case of Eq. ([s]). For filling fraction v — 1/q, q odd, this 
Hamiltonian is given by: 



q-2 N 

E E 

m{odd) i<3 



(6) 



This two-body "hardcore" potential penalizes any state 
where two particles have a relative angular momentum 
smaller than q. The MR Pfaffian wavefunction 



MR 



= Pf 



^ * ^ i<j 

is the exact zero-energy ground state of a three body 
Hamiltonian projecting onto electron triplets instead of 
pairpSl. 

In an ideal, strictly two dimensional system, the 
electron-electron interaction is given by the 2D Coulomb 
potential, V{k) = (e^Z/e)(l/fc). The finite extent of an 
experimental quantum Hall system in the perpendicu- 
lar direction will alter the ideal 2D interaction, yield- 
ing an effective quasi-2D electron-electron interaction. 
There ar e sev eral models for the effect of the finite layer 
thicknes^2ll2lHlll^ how ever, these models all provide sim- 
ilar qualitative result^^^M!^ Therefore, we will focus on 
one particular model, the infinite square well potential. 
In this model, we average the three dimensional Coulomb 
potential over the single-particle ground state n{z) of an 
infinite square well (i.e., n{z) = ^J 2 / d cos{'k z / d)) in the 
perpendicular dimension, yielding the effective interac- 
tion potential 



VsQ{k) 



dzidz2\n{zi)\'^\n{z2)\ 



2 -k\zi-Z2\ 



{^kd+'-^ 



327r"'(l- 



ek 



(kdY[(kdy 
{kdf + 47r2 



(8) 



Combining Eq. ^ with Eq. Q gives us effective pseu- 
dopotentials as a function of finite thickness d. Landau 
level n, and relative angular momentum m. 

We diagonalize the FQH Hamiltonians (one for each 
d/l and LL index n) in the spherical geometrjl^ where 
N electrons are confined to the surface of a sphere. Al- 
though we use this geometry, we use the pseudopotentials 
obtained from the infinite planar geometry (Eq.Q) since 
the finite layer thickness effect is more conveniently mod- 
eled in this case. Furthermore, the pseudopotentials in 
the spherical geometry equal those in the planar geom- 
etry as the thermodynamic limit is approached (as the 
spherical radius is taken to infinity) and it can be argued 
that they provide a better approximation to the thermo- 
dynamic limit (this is discussed in detail in Ref. |43| and 
in the appendix) . In the spherical geometry, the perpen- 
dicular (radial) magnetic field is provided by a magnetic 
monopole of strength, Q, quantized in half-integer units, 
placed at the center of the sphere. The eigenvalues of the 
squared magnitude and z-component of the angu- 
lar momentum, S{S + 1) and m respectively, are good 
quantum numbers for the single particle wavefunctions, 
where S is related to the LL index by the constraint, 
S = \Q\ -\-n and m is constrained such that —S < m < S. 
Thus the degeneracy of a LL with index n is given by 
g = 2S + 1^ 2{\Q\ -I- n) -I- 1. The fiUing factor for a LL 
is defined in the thermodynamic limit v = lim7v->.oo N/g. 
The uniform ground state has total angular momentum, 
L — and therefore, can be obtained in the Hilbert sub- 
space where the total z-componcnt of angular momen- 
tum, Lz = 0. 

In this study, we consider the FQH ground states at the 
Laughlin filling fractions = 1/3 and v — 7/3 with par- 
ticle number N — 6, 7, and 8 and the even-denominator 
filling fractions v = 1/2 and v = 5/2 with particle num- 
ber N = 8 and N — 10. We restrict ourselves to these 
relatively modest system sizes in order to investigate a 
large number of FQH ground states for various values 
of the finite thickness d/l with reasonable computing re- 
sources. Although the Hilbert space for particle number 
iV = 12 at the half fillings is not prohibitively large, this 
system is also ahased with = 2/3 and could, therefore, 
yield ambiguous results. Since we are largely concerned 
with the qualitative features of the finite-thickness effect, 
these system sizes are adequate. 

We calculate the entanglement entropy (EE) and the 
entanglement spectrum (ES) of FQH ground states by 
dividing the sphere into two regions. In particular 
we write our Fock space as the tensor product of two 
subspaces H = Ha <8) 'Hb with respective basis states 

I^a) = 



-S' ' 



e Ua and jt/'l) = 

g_-^B^,.,g_j^B^^^...,,g I e Hb where n„j is the occu- 
pation number of the Landau orbital with angular mo- 



mentum, m and N^, 



iVoib = 25" -hi is the total 
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cial choice of unnormalized basis states. The normaUzed 
single particle wavefu nction on a sphere with angular mo- 
mentum m is given 



f2S + lf 25 



1/2 



where u = cos(6'/2)e*'^/2 ^nd v = sin(6'/2)e"*'^/2 ^j^^ q 
and (f) the usual spherical coordinates. In the conformal 
limit, we "unnormalize" the single particle wavefunctions 
by removing the prefactor in Eq. Q such that the wave- 
functions take the simple form '^{u,v) = 
This procedure is an attempt to remove the finite size 
effects inherent in these calculations by basically remov- 
ing the "length" in the problem. With the ground state 
redefined in this new basis, the entanglement spectrum 
is calculated as described above. 



FIG. 1. Graphical illustration of the partitioning of the Fock 
space. In the spherical geometry, the single-particle states are 
states with the z-component of angular momentum from S to 
— S, represented by the solid latitudinal lines. We choose our 
partitions to cut the sphere in two as close to the equator as 
possible, represented by the dashed lines. Thus, for A'^oib = 
25-1-1 even, we cut the sphere after A'orb /2 (see top panel) 
and after (A^orb -|- l)/2 for A^'orb odd (see bottom panel). 



number of Landau level orbitals. For all cases we choose 
our partitions when dividing our Fock space H into Ha 
and Hb such that for the number of single particle or- 
bitals A^orb ^2S + 1 even, N^^^^ = N^^^^ = N^rh/^ and for 
TVorb odd, N^^^ = A^i?b + 1 = (A^orb + l)/2. Geometrically 
this is equivalent to dividing the sphere along a line of 
latitude (see Fig. [I]). 

We calculate three quantities of interest: 

(i) The bipartite entanglement entropy EE of the 
ground state is given by Se = Tr[p^ log(pA)] where 
Pa is the reduced density matrix defined by pA = 

TTBm in- 

(ii) The entanglement spectrum is obtained from 
the eigenvalues of the reduced density matrix, pi, by the 
simple relation = 21n(|pi|). This is equivalent to find- 



ing the Schmidt-decomposition of the matrix Wij where 



0^) e Ha and 



E.exp(-e./2)|0ii) 



\(f>A), \4'a) ^ '^A and 10^) e He- Since the quantum 
numbers for angular momentum and particle number in 
each region, Lf , , , , are constrained such that 



L, ^ and 



N, the reduced 



density matrix is block diagonal with and being 
good quantum numbers for the eigenstates of pA- There- 
fore Lf and A'^'^ are good labels for the corresponding 
entanglement spectrum. 

(iii) The "conformal limit" of the entanglement spec- 
trum. Recently, Thomale et alW^ have introduced the 
conformal limit when calculating the entanglement spec- 
trum of spherical FQH ground states. This limit is ob- 
tained by expressing the ground state in terms of a spe- 



III. RESULTS 
A. Entanglement Entropy 

We now report numerical results for the entanglement 
entropy (EE) of quasi-2D FQH ground states as a func- 
tion of the finite layer thickness d/l for FQH states in the 
LLL = 1/3 and 1/2) and the SLL {v = 2 + 1/3 = 7/3 
and 2 -f 1/2 = 5/2). As mentioned above we choose the 
partition to be as close to the equator of the sphere as 
possible to minimize finite size effects. 

The results for EE for the Coulomb ground state at fill- 
ing fractions v = 1/3 and v = 7/3 are shown in Fig. [2] as 
a function of finite layer thickness d/l. For comparison, 
the EE of the corresponding Laughlin model wavefunc- 
tion is also shown as a d/l independent horizontal line. 
In each of the figures, we see that in the LLL, the EE is 
near that of the Laughlin model wavefunction at d = 
and rises slightly as a function of d/l. In contrast, the 
EE in the SLL is large compared to that of the Laughlin 
model at d = 0, but decreases as a function of d/l and 
evidently reaches an asymptotic value. The qualitative 
behavior is independent of system size. If we consider 
ASe — Se ~ Se, model for both the LLL and SLL filling 
fractions and speculate that ASe is a qualitative measure 
of how far removed the ground state is from the Laughlin 
model state, then we see that the LL dependence of ASe 
as a function of d/l behaves qualitatively similar to that 
of the overlap between the ground state and the model 
wavefunction as reported in Ref. 1371 IMl SH and |331 In 
particular, the ground state in the LLL is a "strong" 
FQHE state (i.e., ASe is small) at d = and grad- 
ually becomes "weaker" for increasing d/l (albeit only 
slightly), whereas in the SLL, the ground state is ini- 
tially weak at d = but gets stronger with increasing d/l 
(i.e., ASe decreases). Thus, the EE for these cases qual- 
itatively and semi-quantitatively captures how well the 
states are "Laughlin-like" as a function of d/l in similar 
manner to the overlap. 
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FIG. 2. (Color online) Entanglement entropy Se as a function of finite layer thickness, d/l for Laughlin filling fractions u — 1/3 
and v = 7/3 for particle number A*' = 6, 7 and 8. The dashed and dotted lines in the left panels correspond to the Coulomb 
Hamiltonian of a quasi-2D system in the LLL and SLL, respectively, whereas the solid lines in the left panels corresponds to 
the finite size Laughlin states. The plots in the right panels give one minus the percentage difference in the Coulomb EE and 
the model state EE, i.e., 1 — |A5'_e|/S_e, model in the LLL (dash-dotted line) and the SLL (dotted line) and are found to be 
similar qualitatively and quantitatively to overlap calculations^^ ■^^-ii, 



Our operational definition of "weak" and "strong" de- 
pends on how close the EE of the Coulomb state is 
to the model state which, in this case, is the Laugh- 
lin state. In the SLL, Se becomes closer to the Se 
for the Laughlin state but, as mentioned above, appears 
to saturate at some asymptotic value that is still nearly 
^ l-15'£;_Laughiin- In contrast, the EE in the LLL is almost 
identical to that of the Laughlin state. We conjecture as 
to the reason for this difference between the EE in the 
SLL Coulomb ground state compared to the Laughlin 
state and the difference between the EE in the LLL as 
compared to the SLL: (i) it is possible that the FQHE at 
7/3 is not described by the Laughlin state and is instead 
described by a state in a different topological universality 
class such as those given by Read and RezayP^l and Bon- 
derson and SlingerlancPi', (ii) perhaps composite fermion 
interactions, which are thoughlP^ to be more relevant in 
higher LLs, are producing this difference in Se and the 
Laughlin state, (iii) perhaps the 7/3 FQHE state is in 
fact a Laughlin state but our model system is leaving 
out realistic effects such as LL mixing which are crucial 
to its success. 

Fig. [3] gives results for the EE of FQH ground states 



with even denominator filling fractions v — 1/2 and 5/2 
as a function of finite layer thickness d/l. Also shown in 
the figure is the EE of the Moore- Read Pfaffian state for 
comparison. In the LLL {1/ = 1/2), the EE has a weak 
minima as a function of d/l, in contrast to the Laugh- 
lin fractions (this minimum is difficult to discern on our 
scale). The location of this minima changes with N, sug- 
gesting a finite size effect, but the qualitative behavior 
is similar in both cases. In the SLL {v = 5/2), the EE 
has a very pronounced minima that approaches the EE 
of the MR Pfaffian model state for N = 10 and crosses 
it for = 8. This suggest that the FQH states becomes 
more MR Pfaffian-like at near an optimal d/l. However, 
this optimal d/l also changes with N. Similarly to the 
Laughlin fractions, this LL dependence in EE as a func- 
tion of d/l is also qualitatively similar to that seen in the 
overlap between the FQH ground states and MR Pfaffian 
reported in Ref. |42| and |43l These results suggest that 
= 1/2 is not particularly well-described by the MR 
Pfaffian, whereas v = 5/2 is better described by the MR 
Pfaffian model state at finite thickness. 

We note that recently, entanglement entropy in the 
SLL including finite thickness effects has been investi- 
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FIG. 3. (Color online) Entanglement entropy Se as a function of finite layer thickness, d/l for even-denominator filling fractions 
u = 1/2 and u — 5/2 for A'^ = 8 and 10. Similarly to Fig.[2]tlie plots in the right panels give one minus the percentage difference 
between Se and S'b, model- 



gatecPMl], However, the previous study did calculations 
using the torus geometry, in contrast to our spherical 
geometry, and attempted to isolate and calculate the so- 
called topological term of the entanglement entropy. The 
entanglement entropy can be essentially divided into two 
pieces 

SE^aL~-f + o(^^^ (9) 

where L is the linear length of the boundary dividing 
the system into parts A and B (in our case it would 
be circumference of the sphere where we made our cut). 
The aL term is non-topological in origin and the —7 
term is the topological entropy and for the Laughlin 
and MR PfafSan state can be calculated analyticalljEsl 
7 = In -y/m for the Laughlin state and 7 = In ^/4m for the 
MR Pfaffian state at f = 1/m. Extracting the topolog- 
ical entropy from a Coulomb Hamiltonian requires one 
to numerically calculate the exact ground state for many 
different systems sizes and system cuts and perform a 
thermodynamic extrapolation. This is a labor intensive 
procedure that inherently induces statistical errors. Such 
numerical extrapolation, without some strong theoretical 
guidance about the finite-size behavior of the system, is 
often unreliable in estimating quantities in the thermo- 
dynamic limit. 

The conclusion of Refs. [83l and [84l was that the topo- 
logical entropy of the ground states of the LLL or SLL 
Coulomb Hamiltonians was consistent with associated 



model states (we note, however, that in Ref. |83j it was 
concluded that j/ = 7/3 was more consistent with the 
k = A Read-Rezayi statd^ instead of the Laughlin state) . 
However, they also included finite thickness in the form of 
an infinite square well potential and, interestingly, found 
that there was not much difference between the EE and 
the topological entropy with or without finite thickness 
included. We, however, clearly see a finite thickness ef- 
fect on the EE. It is possible that this difference in the 
two studies (our present study and Refs. 1551 and [H?|) is due 
to the different geometry used in the calculations (sphere 
vs. torus) but we find this scenario unlikely since most 
quantities of i nterest produce consistent results in the 
two 

geometriePE^m. Such a comparison between ge- 
ometries (torus vs sphere) was shown in Ref. .86, to give 
similar results for the entanglement spectra of Laughlin 
states, supporting our suspicion. Moreover results given 
in Ref. 60 suggests that the the extrapolation procedure 
performed in Refs. [83] and [84] may have been inappro- 
priate for the torus. More work is clearly necessary to 
understand the difference between the results in spheri- 
cal and toroidal geometries, particularly in the presence 
of the realistic finite thickness effects. 

Before moving on to entanglement spectra we briefly 
discuss how our results compare to the previous overlap 
calculations done in Refs. |42| and l43l The right panels in 
Figs. [2] and [3] gives one minus the percentage error in the 
entanglement entropy, 1 — \^Se\/ Se, model- In Ref. ,43 it 
is found that the overlap between the Laughlin state and 



the Coulomb ground state at 1/3-filling in the LLL and 
SLL is approximately ~ 0.99 a.t d/l = and is reduced 
monotonically to ^ 0.98 at d/l = 8 in the LLL and is 
~ 0.73 at d/l = and has a maximum of ~ 0.84 for d/l ^ 
4 in the SLL. These overlap trends are very consistent 
with what we have seen p reviou sly in EE. For the 1/2- 
fiUed LLL and SLL we fincP^M! the overlap is relatively 
constant in the LLL at ~ 0.9 and in the SLL it is ~ 0.96 
at d/l — and has a maximum value of nearly ~ 1 at 
d/l ~ 4. Again, one minus the percentage error in the 
entanglement entropy tracks the behavior in the overlap 
to a remarkable degree. Perhaps this is not a surprise 
since if the overlap (^ol^modei) is close to one then the 
EE (which is a particular combination of should 
also be nearly identical to the EE of the model state 

^modol- 



B. Entanglement Spectrum 

In the previous section, we saw that the entanglement 
entropy Se (and in particular, A^^;) as a fun ction of d/l 
behaves qualitatively similarly to the overlafP^liS] Yot 
the half-filled case, increasing d/l makes the calculated 
Se closer to the MR Pfaffian state for the SLL {v = 5/2) 
in a rather dramatic way while increasing d/l has very 
little effect on the Se in the LLL = 1/2), i.e., using the 
entanglement entropy as a measure we see that the MR 
Pfaffian is stabilized by finite thickness. For the 1/3-filled 
case we find that increasing d/l drives Se away from the 
Laughlin value in the LLL {i/ — 1/3) and closer to the 
Laughlin value in the SLL (z/ = 7/3), however, as in the 
previous overlap investigations, the value of the entan- 
glement entropy for the 7/3 case never gets as close to 
the Laughlin value as the 5/2 entanglement entropy gets 
to the MR Pfaffian. As discussed above, this could be 
a hint that something is missing from our understanding 
of the physics for the FQHE at v = 7/3. 

To gain a deeper understanding of entanglement, we 
now turn our attention to the finite layer thickness de- 
pendence of the entanglement spectrum (ES), which as 
discussed earlier, provides more information than the EE 
alone. To calculate the ES, we partition the sphere the 
same as was done for the EE. We follow the convention 
established by Li and Haldan^^U and restrict ourselves to 
the part of the ES where the number of particles in the 
A partition, N^, is the same as that of the "root" con- 
figuration for the correspon ding L aughlin or Moore- Read 
Pfaffian model wavefunctiorpl^' for a given partition size 
N^y^. The "root" configurations describe the occupancy 
of LL orbitals for MR Pfaffian and Laughlin model states 
in the thermodynamic limit. Root configurations with a 
maximum z-componcnt of angular momentum, and their 
corresponding quantum numbers, and L^, are given 
in Table[T]for different filling fractions and partition sizes. 

In order to obtain a general qualitative picture of how 
the ES changes as a function of the finite layer thick- 



8 



r lo^ri state 




root conng. 




iV 


Laughlin 1/3,7/3 
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'10010010' 


13.5 


3 




10 


'1001001001' 


18 


4 




11 


'10010010010' 


24 


4 


MR Pfaffian 1/2, 5/2 


7 


'1100110' 


16 


4 




9 


'110011001' 


24.5 


5 



TABLE I. Root configurations of the Laughlin {v = 1/3, 7/3) 
and MR Pfaffian (v = 1/2, 5/2) wavefunctions for the given 
partition sizes, A^^bi on the sphere. '10010010', for example, 
means that the single-particle angular momentum S = 7.5, 
S — 3 = 4.5, and S — 6 = 1.5 are all occupied with the others 
unoccupied. Hence, there are A'^^ = 3 electrons with total 
2-component of angular momentum = 35 — 9 = 13.5 in 
this root configuation. 



ness, we calculate the "entanglement gaps" in each ES 
and plot it as a function of d/l. An entanglement gapl^ 
is defined as the difference between the low-lying "CFT 
levels" (i.e., those levels displaying the low lying CFT 
counting structure) and the higher non-CFT levels for a 
given value of in the spectrum. According to the Li 
and Haldane conjecture, the state has an underlying CFT 
if the entanglement gaps are finite in the thermodynamic 
limit. However only the entanglement gaps at relatively 
"small" values of AL = Lf ,.^^^ — are relevant due to 
finite size effects, where Lf^.^^^ is the total ^-component 
of angular momentum of the root configuration. The fi- 
nite number of LL orbitals limits the number of possible 
"edge excitations." Therefore only a few levels are ex- 
pected to have the same counting structure as the CFT 
edge modes. The "depth" (i.e., the max AL) at which 
the counting structure in the ES is consistent with the 
CFT edge modes is dependent on the system size, N. It 
has been conjectured^^ that these "finite size" effects of 
the entanglement spectra contain information such as the 
generalized statistics of the underlying FQH state, how- 
ever, we will not consider such conjectures in this work. 

We can determine a suitable cutoff for AL by exam- 
ining when the level counting in the ES of the model 
states deviate from the expected counting in the thermo- 
dynamic limit. To illustrate this finite size cutoff, we give 
the ES of the Laughlin state in Fig. |4] For the Laughlin 
state, the multiplicity of CFT levels is given by p(AL) 
where p(m) is the partition function of the integer to. 
The first 7 values of p(to), starting with to = are 1, 1, 
2, 3, 5, 7, and 11. In Fig.g) we see that for = 6 and 7, 
the level counting begins to deviate from p(AL) at AL = 
4, and for A^ = 8, the deviation begins at AL — 5. Thus, 
for our study we will focus on the entanglement gaps for 
AL ~ 0, 1, 2 and 3 for the A^ = 6 and 7 Laughlin sys- 
tems, and for A^ = 8, we also examine the entanglement 
gap at AL = 4. 

We determine the finite size cutoff for the entanglement 
gaps of the half-filled FQH states in a similar manner, 
which we now illustrate. The ES for the MR Pfaffian 
model states are shown in Fig. [5] The counting rules 
for the MR Pfaffian model state depend on where the 
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FIG. 4. Entanglement spectrum of the v = l/?> filled Laughlin model state for particle number A'^ = 6 (top panel), 7 (middle 
panel), and 8 (bottom panel). The finite size cutoff used to examine the entanglements gaps is illustrated by the vertical line. 
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FIG. 5. Entanglement spectrum of the v = 1/2 filled MR Pfaffian model state for particle number N = 8 (top panel) and 
10 (middle panel). The finite size cutoff used to examine the entanglements gaps is illustrated by the vertical line and the 
particular partition of the Fock space is given using the Li and Haldane notatioiP^. 
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partition is made, which correspond to choosing one of 
the three sectors of the corresponding CFT^^. For the 
case of iV = 8, the partition along the equator is equiv- 
alent to the ^"[010] partition in Li and Haldane's nomen- 
clature (i.e., a cut between two unoccupied orbitals in 
the root configuration). The CFT level counting for the 
first 4 Virasoro levels of this partition are 1, 1, 3, and 
5. The counting in the MR Pfaffian ES given in Fig. [5] 
with iV = 8 deviates from this structure at AL = 3. For 
= 10, the partition along the equator corresponds to 
P[l|l] (i.e., a cut between two occupied orbitals of the 
root configuration), which has a CFT level counting of 1, 
2, 4, and 7 for the first 4 Virasoro levels. Examining the 
ES in the figure for N = 10, we see this spectrum also 
deviates from the expected counting at AL = 3. Thus 
for the half-filled FQH states we examine in this study, 
we concern ourselves only with the entanglement gaps up 
to AL = 2. 

The entanglement gaps, which we denote as A^ for 
i = AL, are calculated by finding the difference between 
the largest suspected low-lying CFT level and the next 
highest level at the given value of Lf in the ES of the 
numerically obtained Coulomb ground states for varying 
d/l. We also calculate the minimal gap between CFT and 
non-CFT levels for AL < m, which we denote as Ao_m 
where m is the cutoff described above. The minimal gap 
gives us a qualitative measure of how well separated, 
overall, the CFT levels are from the generic non-CFT 
levels. The suspected low-lying CFT levels are identified 
by the expected counting described above. An exam- 
ple of this procedure is shown in Fig. [6) which shows 
the entanglement gaps in the ES of the v — 1/3 ground 
state for d = and N = 7. Note that throughout this 
work, when presenting figures showing ES, we color code 
the suspected CFT levels with red (gray) diamonds con- 
nected by a red (gray) dash and all other ES levels with 
a black dash. The suspected CFT levels are chosen by 
calculating the ES for the model state (be it the Laugh- 
lin or the MR Pfafhan) and noting how many ES levels 
n{L^) there are for each L^. Then, when we consider the 
ES for the Coulomb Hamiltonians, we identify suspected 
CFT levels (and color code them) as the lowest n{Lf) 
ES levels for each Lf. 

Entanglement gaps as a function of finite layer thick- 
ness d/l for the Laughlin filling fraction = 1/3 arc 
shown in Fig. [7j The entanglement gaps are slightly de- 
creasing with d/l for all cases, indicating that the states 
are weakening. These trends are similar to those ob- 
served in the EE at v = 1/3. Note that the minimal gap 
for = 6 and A'^ = 8 is initially small and becomes zero 
for d/l > 4. This may indicate that the FQH state col- 
lapses at a finite thickness, as has been shown in previous 
works'^* (the previous works showed the FQHE to col- 
lapse at very large d/l). However, this effect is not seen 
in the minimal entanglement gap for N = 7. This "even- 
odd" finite size effect is likely due in part to a trade-off 
between the finite size cutoff and the number of avail- 
able orbitals. Indeed, the finite size cutoff is the same 



for iV = 6 and N = 7, but the larger Hilbert space for 
the N = 7 case allows for more "edge excitations" that 
strengthen each entanglement gap, not just the minimal 
gap, compared to = 6. In all cases, however, the over- 
all trends in the entanglement gaps (i.e., slight decrease 
with d/l) are qualitatively similar to those seen in the 
EE (in particular, ASe) and the overlap in Refs. H^and 

m 

To illustrate this overall trend in the entanglement 
gaps for v = 1/3, we provide the ES of the ground states 
in Fig.jsjfor d/l = 0, 2 and 6. We have marked the levels 
that are consistent with the counting found in the ES of 
the Laughlin model state shown in Fig. |4] for all values 
of AL and indicate our chosen finite size cutoff. We see 
that qualitatively, the ES is largely insensitive to finite 
d/l. Moreover, on the right of the finite size cutoff, except 
for the largest CFT state at AL 4, the low-lying CFT 
levels are well-separated from the higher energy generic 
levels. 

We now examine the case when i/ = 7/3 in compari- 
son. In Fig. [9] are the entanglement gaps as a function 
of finite layer thickness. For A^ = 6, the root entangle- 
ment gap, Ao, generally increases with d/l. Ai and A2 
each have a weak, local maxima near d/l ^ 4 and A3 
is actually decreasing with d/l. Moreover the minimal 
entanglement gap is zero throughout. The entanglement 
gaps for A^ = 7 are are each monotonically increasing 
with d/l, similarly to Aq in the N = 6 case. The min- 
imal gap, which is initially zero, opens at d/l ~ 3 and 
then gradually increases with d/l in this case. The case 
when A^ = 8 shows trends similar to the A^ = 6 case. 
Here Aq, Ai, and A2 increase with d/l, A3 has a local 
maxima near d/l ~ 4, and A4 decreases with d/l. The 
minimal gap for A^ = 8 is zero throughout. Again we 
see an "even-odd" finite size effect in the entanglement 
gaps as was seen with i> — 1/3. However, in this case, we 
have entanglement gaps that increase, decrease, or have 
a weak maxima as a function of d/l. This is in contrast 
to the = 1/3 case where all entanglement gaps follow 
the same trend with finite d/l. The different trends in 
the entanglement gaps may suggest that the topological 
signature of the v — 7 /Z state differs from that of the 
Laughlin state. 

Some illustrative examples of ES at v — 7 /Z are given 
in Fig. [To] with N = % and d/l = 0, 4 and 6. The given 
ES appear to have structure similar to that seen in the 
V = 1/3 case, however, we see that for AL = 4, the 
higher energy "CFT" states are virtually indistinguish- 
able from the "generic states" . This "blending" appears 
to get worse for larger d/l. Again, these results may sug- 
gest that the Laughlin model state is not an accurate 
description for the ly = 7/3 state. 

Results on the entanglement gaps for the even denom- 
inator filling fraction = 1/2 are shown in Fig. 11 Here 
we see that the entanglement gaps slightly decrease with 
d/l and behave similarly to the EE at this filling frac- 
tion. Also note that for N = 8, the minimal gap is small 
and decreases with d/l, while for A^ = 10, the minimal 
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FIG. 6. (Color online) Entanglement spectrum for v — 1/3, d — 0, and A'^ = 7 as a function of ^-component of angular 
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FIG. 9. (Color online) Entanglement Gaps for the Coulomb Hamiltonian as a function of finite layer thickness, d/l for filling 
fraction p = 7 /?> and particle number A'' = 6 (top panel), N — 7 (middle panel), and N — 8 (bottom panel) with partition at 
the equator. 



gap is zero throughout. As mentioned earlier, there has 
been no definitive experimental observation of FQHE at 
^' = 1/2 in monolayer systems consistent with our calcu- 
lations. The full ES of the ground states with v — 1/2 
are given in Fig. 12 for d/l — 0, 2 and 6 and = 8. 
Qualitatively, we see that the ES is largely insensitive 
to the finite-thickness effect. Moreover, the largest sus- 
pected CFT level for AL = 2 is well separated from the 
other CFT levels and appears to be more consistent with 
the generic levels. Again, this suggests that v = 1/2 is 
not described by the MR PfafRan wavefunction. 

Fig. [13] shows the entanglement gaps at filling frac- 
tion 1/ = 5/2. For A^ = 8 each entanglement gap peaks 
at a certain value for d/l. In particular A2 peaks near 
d/l ~ 2.5; the other gaps peak near d/l ~ 4. We also 
see peaks in the entanglement gaps for the case when 
A^ — 10. Here, the gaps gradually rise to a local maxima 
near d/l ~ 1.5 and then slowly decay for increasing d/l. 
Note that the gaps in this case are generally smaller com- 
pared to those observed for A^ = 8. These results may 
suggest that there is a slight difference in the finite-size ef- 
fect on the different MR PfafRan CFT sectors. However, 
these results are qualitatively similar to the EE results 
and the results on the overlap in Refs. |33|(i.e., the MR 
Pfaffian signature of the v — h/2 state is strengthened 
by the finite size effect). 

We also provide the ES of the v — b/2 state for A = 8 



in Fig. 14 for d/l — 0, 4, and 6. Here, we see the ES 
"opens" at d/l — 4, giving a larger separation between 
the CFT and generic levels in the spectrum compared to 
d/l = and 6. Again, these results suggests the ly — 5/2 
is, indeed, described by the MR Pfaffian wavefunction, 
and this description is more stable at finite thickness. 

In summary, the entanglement gaps in the ES have a 
similar dependence on finite thickness as the EE, leading 
to similar conclusions. However finite size effects pre- 
vent us from making definitive statements. In the next 
section, we attempt to alleviate this problem using the 
conformal limit. 



C. The Conformal Limit 

In the previous section, we used entanglement gaps in 
the ES to evaluate the "strength" of a state as a function 
of d/l and we were able to confirm the MR PfafRan sig- 
nature of the I' — 5/2 state and distinguish it from the 
(lack of) signature of the 1/ = 1/2 state. However, we 
have intentionally ignored a significant part of the ES in 
order to avoid finite size effects, i.e., we focused on the 
region of the ES with small AL (see Fig. [6|. We deter- 
mined the size of this region by examining where the ES 
of the finite sized MR Pfaffian and Laughlin model states 
deviate from the conjectured structure in the thermody- 
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FIG. 10. (Color online) Entanglement spectrum for the Coulomb Hamiltonian as a function of z-component of angular 
momentum for filling fraction v = 1/3 and particle number N = 8 for d/l = (left panel), d/l — 4 (middle panel), d/l — 6 
(right panel). The suspected CFT levels consistent with the Laughlin model state for each are marked. 



namic limit (i.e., the edge state level counting given by 
the suspected CFT). Thus, we have confirmed the MR 
PfafRan and Laughlin signatures only to a certain extent 
because, in fact, there is not an actual entanglement gap 
in the ES. 

It has been conjectured that the full entanglement 
spectrum of the finite sized model states contain informa- 
tion on the topological signature of the FQH state.'^^Ml 
Thus, all states in the ES can be used to identify the 
topological quantities. With this in mind, we now exam- 
ine the entanglement spectrum of quasi-2D FQH states 
in the "conformal limit" (CL), which reportedly allows 
the use of the entire spectrum to examine the state by 
unambiguously defining a full entanglement gap. As dis- 
cussed briefly above, and at length by Thomale et al. in 
Ref. [Ml the CL works by removing finite size effects due 
to the curvature of the sphere and gives an ES with a 
"full" unambiguous entanglement gap in the spectrum 
for topologically ordered states. Thus the presence of an 
entanglement gap in the conformal limit is conjectured 
to be a sign of topological order. A demonstration of 
an ES before and after the CL is given in Fig. [15] Af- 
ter taking the CL of an ES (CLES), we determine the 
"minimal gap" by taking the difference between the high- 
est suspected CFT level and the lowest generic level in 
the entire spectrum. The suspected CFT levels are de- 



termined by comparing the CLES to that of the model 
state with the assumption that all levels in the ES of the 
model state are CFT levels. For comparison, we examine 
the entanglement gaps for each value of Lf and define 
the "average gap" as the average of the individual en- 
tanglement gaps. We also define the "maximum gap" as 
the maximum of the entanglement gaps. Individual gaps 
that are near infinite (i.e., no levels above the highest 
CFT level) are ignored. The minimal gap, the average 
gap, and the maximum gap are calculated for each CLES 
as a function of the finite layer thickness, d/l. 

CLES entanglement gaps as a function of finite layer 
thickness d/l for v = 1/3 (LLL) are shown in Fig. 16 For 
A'^ = 6 and 7, the entanglement gap measures decreases 
with d/l but remains finite throughout. This behavior is 
qualitatively similar to the ES gaps for small AL, as well 
as the EE results, suggesting a weakening of the Laugh- 
lin state. The fact that the minimal entanglement gap 
in the ES for the A^ = 6 case (Fig. [t]) differs from the 
minimal gap in the CLES may indicate that the closing 
of the gap in the ES is due to finite size effects related 
to the curvature of the geometry rather than the limited 
number of LL orbitals. However, the minimal gap for the 
case where N = 8 seems anomalous. Although the aver- 
age and maximum gaps follow similar qualitative trends, 
the minimal gap is at or near for all values of d/l, in- 
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FIG. 11. (Color online) Entanglement Gaps for the Coulomb Hamiltonian as a function of finite layer thickness, djl for filling 
fraction v and particle number N — 8 (top panel) and A'^ = 10 (bottom panel) with partition at the equator. 
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eluding d = 0. How to interpret this result is unclear 
since there is a general consensus that the Laughlin state 
does, indeed, model the = 1/3 state. We can shed 
some light on this anomaly by examining the CLES of 



the FQH states directly. Fig. [T7| shows the CLES for the 
ly = 1/3, N — 8 FQH state at finite thicknesses d/l — 0, 



2 and 6. The suspected CFT levels are marked in each 
plot. 

We examine CLES in the SLL case {i^ = 7/3) in Fig.[l8 
In general, the behavior of each gap measure differs with 
varying d/l. The minimal gap appears fragile and virtu- 
ally disappears for larger N. The average gap has two 
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FIG. 13. (Color online) Entanglement Gaps for the Coulomb Hamiltonian as a function of finite layer thickness, d/l for filling 
fraction v = and particle number N — 8 (top panel) and A'^ = 10 (bottom panel) with partition at the equator. 



local maxima in d/l for N = 6. Only one of the local 
maxima in the average gap is preserved when we look at 
the N — 7 case, and for TV = 8, the average gap fluc- 
tuates. The maximum gap, in general, increases with 
increasing d/l but has a notable peak near d/l ^ 0.7 for 
= 8. The inconsistency in these results may suggest, 
from the ES and EE results, that the Laughlin model 
state is not a suitable model for the v = 7/3 state, or 
other ignored effects are needed for the Coulomb state to 
be adequately described by the Laughlin state. 

Fig. [19] shows the CLES resuhs for = 7/3, iV = 8 
at finite thickness d/l — 0, 4, and 6. We note that for 
each value of d/l, there is very little separation between 
the suspected CFT levels and the generic levels. Indeed, 
if the suspected CFT levels were not marked, there is no 
clear entanglement gap across the whole spectrum. How- 
ever, there does appear to be structure in the CLES for 
small values of AL (i.e., near the "root" configuration). 
What this may imply about the topological signature of 
the V = 7/3 state is not clear. 

Results for the CLES gap measures in the even denom- 
1/2 FQH state are shown in Fig. 



20 



We see 



mator v 

that the minimal gap is nonzero and gradually increases 
with d/l for iV = 8. However, for = 10 the minimal 
gap is zero throughout. The maximum and average gaps 
decrease with d/l for N — 8. For TV = 10, the average 
gap has several local maxima, while the maximum gap 
decreases then suddenly becomes constant with d/l with 



two sharp peaks. Given our results on the EE and the 
ES for this state, the inconsistency between the = 8 
and iV = 10 in the CLES gap measures may suggest that 
the MR PfafRan model state is not a suitable model for 
1/ = 1/2. We also provide the CLES of the ground states 
in Fig. [2l] for = 10 and d/l = 0, 4, and 6. Qualita- 
tively, the CLES do not change very much as a function 
d/l, and there is no clear separation between the CFT 
and generic levels. This, again, suggests that there is no 
FQH state at this filling fraction. 

The CLES gap measures for = 5/2 in the SLL are 
given in Fig. [22j For A^ = 8, the minimal gap has a very 
pronounced peak near d/l w 4. The average and max- 
imum gaps, however, have local minima near where the 
minimal gap is maximum. These "cusps" are a result of 
level crossings. For N = 10, the minimal gap is initially 
zero, but becomes finite for non-zero d/l and peaks near 
d/l « 4.5. The average and maximum gap in this case 
have similar shapes with a peak near d/l « 3. These 
results are qualitatively similar to the results of the ES, 
EE, and the overlap in Refs. EH and |33] Moreover, the 
difference between N = 8 and A^ = 10 may suggest that 
finite thickness affects each sector of the suspected CFT 
differently, but larger system sizes are necessary to verify 
this. In the CLES plots shown in Fig. [23] for A^ = 10 



and d/l = 0, 4 and 6 respectively, we see the entangle- 
ment gap between CFT and generic levels "open" at finite 
d/l = A compared to d/l = and 6. These results are 
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FIG. 14. (Color online) Entanglement spectrum for the Coulomb Hamiltonian as a function of z-component of angular of 
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FIG. 16. (Color online) Conformal limit entanglement gaps for the Coulomb Hamiltonian as a function of finite layer thickness, 
d/l for filling fraction f = 1/3 and particle number A'' = 6 (top panel), N = 7 (middle panel), and N = 8 (bottom panel). 
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FIG. 17. (Color online) Conformal limit entanglement spectrum for the Coulomb Hamiltonian as a function of z-component of 
angular momentum Lf for filling fraction v = 1/3 and particle number N = 8 for d/l = (left panel), d/l — 2 (middle panel), 
d/l — 6 (right panel). The suspected CFT levels consistent with the Laughlin model state for each are marked. 



consistent with results observed with the EE and the ES, 
indicating that MR Pfaffian signature strengthens with 
a finite d/l. 

In summary, taking the conformal limit of the entan- 
glement spectra provides us with a full entanglement gap 
in most cases with a finite thickness dependence that is 
qualitatively similar to the results on the EE. The no- 
table exceptions are the ly = 1/2 which has little or no 
entanglement gap consistent with experimental observa- 
tions, the = 1/3 case at = 8 that is not expected 
given results with other system sizes, and the v — 7/3 
case that is consistent with the results on the EE and 
ES suggesting that other physics besides the Laughlin 
state alone is needed to explain this FQHE. The case 
with N = 8 and v = 1/3, however, is inconsistent with 
most theory and experiment, but when we examine the 
spectra directly, there are a few "spurious" states that 
cross an otherwise full gap. The origin of these "spu- 
rious" states are related to our use of planar Haldane 
pseudopotentials rather than spherical pseudopotentials 
and is discussed in the appendix. However we do not 
expect this choice to alter the topological features of the 
state. Therefore, this result may suggest that a "full" 
quantitative entanglement gap is not necessary to iden- 
tify a topological state. In the next section, we introduce 
the concept of a entanglement spectral density of states 



where a qualitative, "soft" gap may be identified in such 
cases. 



D. Entanglement Spectrum Density of States 

In the entanglement results presented above, we re- 
quire a model state wavefunction for comparison in order 
to systematically define and calculate the entanglement 
gaps. These methods have the obvious disadvantage of 
requiring an ansatz for comparison. In the conformal 
limit case, we assume the low-lying states in the entan- 
glement spectrum should have the exact counting as seen 
in the model entanglement spectrum. This assumption 
may be premature since other finite size effects may cause 
the counting to deviate, even after taking the conformal 
limit especially at the largest AL, see Fig. [TTj With this 
in mind, we attempt to obtain a general qualitative sense 
for how the entanglement spectra vary with finite layer 
thickness by extending the analogy with "energy levels" 
a bit further by calculating the "density of entanglement 
spectral states." With the density of states, we can qual- 
itatively look for entanglement gaps without relying on 
a model state for comparison. Also, we may be able to 
detect "soft" gaps where a small number of states may be 
present within an otherwise prominent gap between two 
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peaks in the density of states. Thus in this section we 
briefly examine this extension by providing results for 
the density of states (DOS) of the entanglement spec- 
trum, both with and without the conformal limit, as a 
hmction of finite layer thickness, d/l. 

The plots shown in Fig. [24] give the density of states of 
the ES and CLES for = 1 /3 as a function of finite layer 
thickness d/l for N = 7. In the DOS for the ES before 
taking the CL, we see sparse low lying states that are sep- 
arated from a denser cloud of higher states by a series of 
gaps. These low lying states are the CFT states from the 
Li and Haldane conjecture. The states appear, largely to 
be insensitive to the finite layer thickness. Turning to the 
DOS for the CLES, a clear gap is much more evident for 
the N = 7 cases. Here the higher-energy states appear to 
makeup a wide, low-density band that is well-separated 
from a low, dense band of states by a gap that decreases 
with d/l. This case seems to illustrate the effect of using 
the CL. 

For comparison, we provide the DOS results for = 8 
in Fig. 



25 



Here there also appears to be low-lying CFT 
states in the ES below a high density region of higher 
energy states. In the DOS of the CLES, a "clear" gap 
does not appear. But the low lying band in this case does 
appear qualitatively similar to the N = 7 case. One may 
possibly associate a "soft" gap in this case, where a few 
states appear to be present between two somewhat dis- 
tinct regions in the DOS. This "soft" gap is qualitatively 



similar to the "clear" gap in the N — 7 case and it does 
appear to slightly decrease as a function of d/l. However, 
it is difficult to distinguish this "soft" gap from the other 
small gaps in the spectrum. 

In Fig. [26] we provide DOS plots for the z/ = 5/2 FQH 
state for iV = 10 as a function of d/l. In this case, the 
ES is especially sensitive to finite layer thickness. How- 
ever, we still see a series of small gaps separating thin, 
dense bands at lower energies. After taking the CL, a 
clear gap at finite (non-zero) thickness has a definite peak 
corresponding to a level crossing. Below the gap, there 
appears to be some band crossings as d/l is varied. 

In summary, the DOS of the entanglement spectra 
(with and without the conformal limit) gives us a general 
qualitative picture of how the ES evolve with a varying 
parameter (i.e., the finite layer thickness d in our case). 
Thus we expect the DOS of the ES to be a good initial 
cursory tool in examining topological states with varying 
parameters. 



IV. CONCLUSIONS 

In this work we study entanglement in finite sized, 
quasi-2D FQH states via the entanglement entropy and 
the entanglement spectrum as a function of the finite 
layer thickness of the transverse dimension in a realistic 
FQH system and compare them to the entanglement sig- 
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FIG. 19. (Color online) Conformal Limit Entanglement spectrum for the Coulomb Hamiltonian as a function of z-component 
of angular momentum Lf for filling fraction u = 1/3 and particle number N = 8 for d/l — (left panel), d/l = 4 (middle 
panel), d/l — 6 (right panel). The suspected CFT levels consistent with the Laughlin model state for each are marked. 
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FIG. 20. (Color online) Conformal Limit Entanglement Gaps as a function of finite layer thickness, d/l for filling fraction 
1/ = 1/2 and particle number N — 8 (top panel) and A*' = 10 (bottom panel) 
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FIG. 21. (Color online) Entanglement spectrum as a function of z-component of angular momentum for filling fraction 
1/ — 1/2 and particle number N = 10 for d/l = (left panel), d/l = 4 (middle panel), d/i = 6 (right panel). The suspected 
CFT levels consistent with the MR Pfafhan model state for each are marked. 




FIG. 22. (Color online) Entanglement Gaps as a function of finite layer thickness, d/l for filling fraction u — 5/2 and particle 
number N = 8 (top panel) and A'' = 10 (bottom panel) with partition at the equator. 
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FIG. 23. (Color online) Entanglement spectrum as a function of z-component of angular momentum for filling fraction 
1/ — 5/2 and particle number A'' = 10 for d/l — (left panel), d/l = 2 (middle panel), d/l — 6 (right panel). The suspected 
CFT levels consistent with the MR Pfafhan model state for each are marked. 
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FIG. 24. (Color online) Density of entanglement energies 
before (top panel) and after (bottom panel) the conformal 
limit for u — 1/3 and particle number = 7 as a function of 
finite layer thickness d/l. 



natures of the Laughlin and MR Pfaffian model states. 
For the Laughlin filling fractions, we find that the EE in- 
creases (decreases) with finite layer thickness for v — 1/3 
{v — 7/3) in the LLL (SLL) with increasing (decreas- 
ing) deviation from the EE of the Laughlin model state. 
However the EE in the SLL reaches an asymptotic value 
larger than the EE of the Laughlin state, possibly sug- 
gesting the v = 7/3 state is modeled by different physics 
than the Laughlin state. Similar behavior is also seen 
in the entanglement gaps of the ES for the Laughlin fill- 
ing fractions. Here we find that the entanglement gaps 
decrease with finite layer thickness for the Laughlin fill- 
ing fractions in the LLL. But in the SLL, the behav- 
ior of the entanglement gaps depend on the "depth" of 
the gap. These results suggest that the Laughlin FQH 
states "weaken" with increasing thickness in the LLL, 
whic h is co nsistent with previous work on quasi-2D FQH 
state^^, but in the SLL, other physics beyond just the 
Laughlin state alone is needed to describe the FQH state. 
The LL dependence of the finite thickness effect at half- 
filhng differs sfightly. The EE of the v = 1/2 state in 
the LLL is largely insensitive to the finite layer thickness 
in contrast to that of the SLL v — h/2 state where the 
EE has a local minima that approaches the EE of the 
MR Pfaffian at finite d/l. This qualitative behavior is 
also seen in the entanglement gaps of the ES for half- 
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FIG. 25. (Color online) Density of entanglement energies 
before (top panel) and after (bottom panel) the conformal 
limit for u = l/Z and particle number A*' = 8 as a function of 
finite layer thickness d/l. 
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FIG. 26. (Color online) Density of entanglement energies 
before (top panel) and after (bottom panel) the conformal 
limit for z/ = 1/3 and particle number A*' = 7 as a function of 
finite layer thickness d/l. 



filled LLs. For = 5/2 in the SLL, we see peaks (local 
maxima) in the entanglement gaps at finite thickness, 
suggesting the z/ = 5/2 is more "MR Pfaffian-like" at an 
optimal thick ness, which, again is consistent with previ- 
ous wory^^EHand strongly suggests the v — h/2 state is, 
indeed, MR Pfaffian. In contrast the entanglement gaps 
of the ES for the v = 1/2 state suggest that it is not mod- 
eled by the MR Pfaffian. Thus, the entanglement gaps in 
the ES allows us to differentiate the v = 1/2 and v = h/2 
states, which we could not definitively establish with the 
EE or the overlap calculations. Of course, we must be 
cautious with these results since the calculated entangle- 
ment gaps made use of only a few "Virasoro levels" in 
the low-lying CFT due to finite-size effects. Assuming 
the Li and Haldane conjecture to be correct, we can only 
say that we have observed the Laughlin and MR Pfaffian 
signatures up to a few "Virasoro levels" . 

We also investigate the conformal limit of the entan- 
glement spectrum which is conjectured to remove curva- 
ture in the spectrum due to finite size effects and allow 
the use of the entire spectrum to determine the topolog- 
ical signature of the state. Our results on the conformal 
limit, however are inconsistent between varying system 
sizes and are difficult to interpret. This appears to be 
due to our choice of using planar pseudopotentials rather 
than spherical pseudopotentials in obtaining the FQHE 
ground states. In the appendix we examine this choice 
by comparing the entanglement spectra of ground states 
obtained by using either spherical or planar pseudopoten- 
tials at d = and observe that the conformal limit can 
be affected by components of the ground state that have 
exponentially small contributions and, therefore, are sen- 
sitive to minor details in the interaction (such as the dif- 
ference between planar and spherical pseudopotentials). 
Thus the presence of the entanglement gap in the con- 
formal limit is sensitive to certain details in the effec- 
tive interaction that may not be relevant in determining 
the topological features of the state. Further work using 
much larger system sizes would be necessary to resolve 
this issue which is well beyond the scope of the current 
work. 

We have also introduced the notion of entanglement 
density of states as a method for examining the idea of an 
entanglement gap without an explicit reliance on a model 
wavefunction. Although, far from definitive, the entan- 
glement DOS suggests itself as a powerful tool to deter- 
mine the topological nature of a particular ground state. 
Our detailed numerical study establishes the entangle- 
ment DOS to be a useful quantity underlying topological 
FQHE particularly in the context of finite size numerical 
calculations. 

It is interesting to observe that the entanglement mea- 
sures give results similar to those obtained with over- 
laps in Refs. and H51 Whereas the overlap is a simple 
measure of how well a numerically obtained ground state 
matches a particular model state (e.g. the Laughlin state 
or the MR Pfaffian state), the entanglement measures (in 
particular, the ES) is a more general measure of how well 
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a state fits a suspected CFT (i.e. universality class) that 
describes the model state. Therefore, it can be said that 
these results confirm the conclusions in Refs. |311 and 03] 
in a more general sense in respect to the Laughlin and 
MR PfafHan CFTs. However, we must be cautious in 
this generalization given that we have only observed the 
Laughlin and MR Pfafhan signature up to a few "Vira- 
soro levels" and different theories can result in the same 
low-level structure in the E^^H More work is necessary 
to understand how well entanglement measures can def- 
initely identify universality classes in finite systems. 

In interpreting our results and conclusions, one may 
wonder about the importance of finite size effects on our 
numerical diagonalization. The possible limitations asso- 
ciated with finite-size effects are of course always present 
in any exact diagonalization study of any FQHE system, 
and the possibility that some of the conclusions are af- 
fected by finite size effects can never be ruled out even if 
the calculations are carried out on systems much larger 
than what we use in this work, since in the end any state- 
ment about an experimental system based on calculations 
performed on few-particle systems is always subject to an 
extrapolation to the thermodynamic limit. We believe 
that all our conclusions regarding the importance of fi- 
nite quasi-2D thickness effect on the FQHE entanglement 
spectra are valid independent of the rather modest size 
of o ur fini te system diagonalization study because earlier 
worlP2l^clearly established, when compared with calcu- 
lations I22IIII carried out on much larger systems, that 
the system size we use in this work, namely N=8, is 
certainly adequate in making qualitatively correct con- 
clusions about the SLL FQHE. Our goal in this paper 
has been to study as many FQHE states as feasible as 
a function of the quasi-2D layer thickness in depth, thus 
necessarily (due to the computational time restrictions) 
limiting our system size to N=8 which should be ade- 
quate. Nevertheless, we feel that future work should ex- 
plore larger system size diagonalization in order to study 
the finite-thickness effect on the entanglement spectra of 
various FQHE states. 

We also note that there are various alternatives to the 
infinite square well effective potential in examining the 
finite thickness effect. The Zhang-Das Sarma potential 
is likely the most well-known and oft-used alternative in 
this regard. We choose to focus on the infinite square 
well instead of the Zhang-Das Sarma potential because 
the infinite square well is more closely aligned to realistic 
quasi-2D systems. Moreover many observables calculated 
with the Zhang-Das Sarma potential besides the overlap 
have been shown to give similar results compared to the 
infinite square well. Therefore we expect the two poten- 
tials to give similar results with entanglement measures 
as well, but since investigating the nature of these poten- 
tials is not a goal of the present study, we leave explicit 
verification of this assertion for a future work. 

It should be noted that the concept of entanglement 
spectra (or for that matter, the entanglement entropy it- 
self) has no direct experimental or observational implica- 



tions since it cannot be directly measured. The concept 
is useful conceptually and theoretically in ascertaining 
the quantum topological nature of a particular interact- 
ing Hamiltonian, and in that sense its experimental con- 
sequences are indirect since the topological nature of a 
system has obvious experimental consequences. We have 
investigated in this work the utility of the concept of 
the entanglement spectra in ascertaining the underlying 
topological nature of realistic FQHE states as a function 
of the quasi-2D layer thickness, finding that the entangle- 
ment spectrum provides results consistent with what has 
earlier been established in the literature based on the 
wavefunction overlap studies. Our very detailed study 
also indicates that at this stage of theoretical develop- 
ment, the entanglement approach is perhaps no more pre- 
dictive in providing experimental implications of various 
FQHE states than what is already available in the liter- 
ature based on the direct wavefunction overlap studies. 
Further work would be necessary to see if the entangle- 
ment approach has some specific advantages in predict- 
ing experimental properties of FQHE states not already 
apparent in wavefunction-based analyses. 

In conclusion, we have extended the concept of topo- 
logical entanglement spectra and entanglement gaps to 
finite-thickness FQH systems by calculating the FQHE 
topological properties systemically as a function of fi- 
nite thickness of the quasi-2D systems, establishing in 
the process that the FQHE entanglement measures cal- 
culated as a function of system thickness are completely 
consistent with the results obtained earlier in the lit- 
erature using wavefunction overlap calculations. While 
our work establishes various entanglement measures as 
important theoretical quantities classifying FQHE, more 
work will be necessary to understand the finite size as- 
pects of entanglement spectra and entanglement gaps in 
the context of realistic fractional quantum Hall systems. 
Although it is gratifying that the qualitative conclusions 
of our entanglement-measure-based results in this work 
are completely consistent with earlier FQHE results ob- 
tained on the basis of wavefunction overlap calculations, 
it remains to be seen whether the entanglement-measure 
based probes have more predictive power regarding the 
nature of FQHE than the wavefunction-overlap based 
probes or it is simply a deeper way of looking at the 
same physics with no obvious additional implications for 
the experimental occurrence of FQHE. 



Appendix A: Planar vs spherical pseudopotentials at 

d = 

The analysis presented above is based on ground state 
wavefunctions obtained by diagonalizing the quasi-2D 
Coulomb potential in a spherical geometry. However the 
Haldane pseudopotentials used to construct the Hamil- 
tonian are derived from a infinite planar geometry rather 
than a spherical geometry. We choose to use planar 
rather than spherical psuedopotentials because (i) the 
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are given in fig. 27 In the figure, we see that the ES 



27 a) is qualitatively 



TABLE II. Overlap integrals between 1) the exact ground 
state wavefunction using spherical (I'I'sphcre)) and planar Hal- 
dane pseudopotentials (l^'piane)), and 2) the overlap between 
the Laughlin or Pfaffian wavefunction (I^I/modoi) and the ex- 
act ground state wavefunction using spherical or planar pseu- 
dopotentials. 



effective Coulomb potential in a quasi-2D system is more 
naturally obtained in the infinite planar geometry and (ii) 
we expect the spherical and planar pseudopotentials to be 
indistinguishable in the thermodynamic limit. Moreover 
given the mostly qualitative nature involved in studying 
entanglement spectra, we expected this choice to make 
little difference in the results. Nevertheless, there are 
cases under study where this choice matters. The goal 
of this appendix is to highlight some of these cases. We 
show that for d = 0, the low energy spectrum in the 
entanglement spectra are qualitatively similar between 
ground states obtained from either spherical or planar 
pseudopotentials, but higher energy spectra can can dif- 
fer in some cases. This difference does not change the 
qualitative conclusions drawn from the low energy spec- 
tra, but when we consider the conformal limit which looks 
for a full entanglement gap, the difference can lead to dif- 
ferent conclusions (in particular, the case of = 1/3 with 
A'' = 8). We leave the comparison of cases with d > 
and larger TV for future work. 

In Table [TI| we provide several overlap calculations be- 
tween exact ground states at d = obtained using ei- 
ther spherical or planar pseudopotentials. In column 3 
of the table, we see that the overlap between the ground 
states from the spherical and planar cases is generally 
high. The notable exception is the case when = 10 
and V — 1/2 where the overlap is 0. In this case the 
ground state obtained with the planar pseudopotentials 
possesses a different symmetry compared to the ground 
state of the spherical case, which leads to a vanishing 
overlap. Excluding these, columns 4 and 5 of the table 
show that the overlap between the spherical and planar 
ground states with the model Laughlin or MR Pfaffian 
states are qualitatively similar. 

We now turn our attention to the entanglement spec- 
tra and how they may differ with choice of pseudopo- 
tentials. ES (with and without the conformal limit) for 
the exact ground state of the FQHE state at = 7 
and = 1/3 using spherical and planar pseudopotentials 



with planar pseudopotentials (Fig. 
similar to the spectra obtained with spherical pseudopo- 
tentials (Fig. 27 b). The same can also be said with the 
ES in the conformal limit between the planar case (Fig. 
27 c) and the spherical case (Fig. 27 d). Thus, given 
the results in Fig. [27j we would expect that the choice 
of pseudopotenials makes little difference in obtaining a 
qualitative understanding of the ES in this case. 

Fig. [28] compares the ES of the FQHE state at fill- 
ing fraction v = l/i with = 8. In this case we see 
that in the ES before the conformal limit (Fig. 28 a 
and b), the low energy spectra are qualitatively simi- 
lar between the planar and spherical cases. The higher 
energy spectra in the ES, however, show notable differ- 
ences with the planar case having a few CFT levels at 
much higher energy compared to the spherical case. In 
the conformal limit, these higher energy CFT levels lead 
to a vanishing entanglement gap in the conformal limit 
for the planar case (Fig. [28] c) compared to the spherical 
case (Fig. 28 d) where there is a full entanglement gap. 
These are the same "spurious" levels identified earlier 
These results suggest that the vanish- 



in section III C 



ing minimal gap seen in Fig. 16 is due to our choice of 
planar rather than spherical Haldane pseudopotentials. 
This may seem surprising given the large overlaps seen 
in Table [ll] However the states associated with the the 
higher energy CFT levels have exponentially small con- 
tributions to the ground state wavefunction, and thus 
contribute little to the overlap. We might also expect 
these states to be more sensitive to certain quantitative 
details of the potential that do not affect the qualitative 
picture of the FQHE ground state (e.g. values of Vm for 
"large" m). Thus when taking the conformal limit, the 
choice of pseudopotential may matter in some cases in 
order to observe a full entanglement gap. But a qualita- 
tive understanding can still be gleaned from the planar 
case since there does appear to be two distinct regions in 
the CLES that we can identify, at least qualitatively, as 
CFT and non-CFT levels. 

We now compare the spherical and planar pseudopo- 
tentials in the SLL with = 7/3. Fig. l29] gives the ES ([29 
a and b) and CLES ([29] c and d) for the v = 7/3 FQHE 
state obtained with either planar or spherical pseudopo- 
tentials with A^ = 8. The planar and spherical cases are 
qualitatively similar in both the ES and CLES and both 
suggest that the Laughlin wavefunction may not describe 
this state, as discussed in sections ]III A[ [1116] and jlll C[ 
Results for the even denominator filling fraction v = 
1/2 with A^ = 10 are given in Fig. 30 In this case, the 
planar results (28 a and c) differ considerably from that 
of the spherical case ( [28] b and d) . This is not surprising 
since the overlap between these two states given in table 
\n\ vanishes. However, it appears that neither state is 
consistent with the MR Pfaffian. 

Comparison of FQHE ground states obtained with pla- 
nar and spherical pseudopotentials for the v = 5/2 state 

1/3 



with A" = 10 is given in Fig. 31 Similar to the v 
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FIG. 27. (Color online) Entanglement spectra and conformal limit entanglement spectra of the exact FQHE ground state for 
= 7 at filling fraction v — 1/3 obtained with either planar or spherical Haldane pseudopotentials at d = 0. CFT states 
associated with the Laughlin model wavefunction are marked. 



case, the low energy spectra in the ES (31 a and b) are 
qualitatively similar between the two cases. The higher 
energy levels in the spectra do differ, but the CLES (31 c 



and d) does appear to give the same qualitative picture. 
Recall that in section [III C[ the minimal gap for this case 
becomes non-zero only at finite d for the planar case. We 
would expect a similar result to occur using the spherical 
psuedopotentials. Verification of this is left for a future 



work. 
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FIG. 30. (Color online) Entanglement spectra and conformal limit entanglement spectra of the exact FQHE ground state for 
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FIG. 31. (Color online) Entanglement spectra and conformal limit entanglement spectra of the exact FQHE ground state for 
= 10 at filling fraction v — 5/2 obtained with either planar or spherical Haldane pseudopotentials at d = 0. CFT states 
associated with the MR PfafBan model wavefunction are marked. 



